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The equilibrium shapes and energies of coherent strained InP islands grown on GaP have been in- 
vestigated with a hybrid approach that has been previously applied to In As islands on GaAs. This 
combines calculations of the surface energies by density functional theory and the bulk deformation 
energies by continuum elasticity theory. The calculated equilibrium shapes for different chemical 
environments exhibit the {101}, {111}, {111} facets and a (001) top surface. They compare quite 
well with recent atomic-force microscopy data. Thus in the InP/GalnP-system a considerable equi- 
libration of the individual islands with respect to their shapes can be achieved. We discuss the 
implications of our results for the Ostwald ripening of the coherent InP islands. In addition we 
compare strain fields in uncapped and capped islands. 
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I. INTRODUCTION 

It was observed about ten years ago that during the 
hetero-epitaxial growth of lattice-mismatched semicon- 
ductors small dislocation-free islands can form.ETO To 
emphasize the zero-dimensional character of their elec- 
tronic density of states, these objects have been labeled 
quantum dots. Since then, the quantum dots have at- 
tracted an enormous interest in the area of semiconduc- 
tor physics.B Besides being fascinating objects for basic 
research, there are also various potential applications, 
ranging from improved device properties of quantum-dot 
semiconductoijJasers, which have in fact already been 
demonstrated!] up to more exotic applications as part 
of a quantum computer.El The "self-organized" hetero- 
epitaxial growth of arrays of quantum dots with a prefer- 
entially narrow size distribution is thus the aim of intense 
research. To finally be able to optimize the growth pa- 
rameters, a detailed knowledge of the growth mechanism, 
its energetics and its kinetics, is essential. 

There is general agreement about the main driving 
force for, e.g., InAs/GaAs or InP/GalnP island for- 
mation in the Stranski-Krastanov growth-mode. The 
growth begins with the deposition of a highly strained 
two-dimensional film, the wetting layer. With the ad- 
dition of more material beyond a critical thickness the 
film becomes meta-stable. Coherent, i.e., dislocation- 
free islands form due to the energy gained by strain 
relaxation in the islands. However, the details of the 
growth mechanismB are not yet well understood. In fact, 
there are competing theories, which are based on energy 
ground-state considerations on the one hand, jOr xtn ki- 
netic, i.e. non-equilibrium effects on the other.EETtfl Ac- 
cordingly, the final destiny of the islands, as predicted by 



the various theories, is rather different. The more con- 
ventionaL&te of the islands would be to undergo Ostwald 
ripening, t3'll3 i.e., the larger islands would grow at the ex- 
pense of the smaller ones. The resulting island size dis- 
tribution would be rather broad. In fact, for the systems 
of concern in our study, Ostwald ripening takes place at 
a slower rate than the growth rate of the islands. Ki- 
netic effects could be active that effectively decrease the 
growth rate aLLhcilarger islands, thus sharpening the size 
distribution.Miljlll-ta--contrast to these growth scenar- 
ios, Shchukin et alx3^ have suggested that there exists 
a range of parameters for which islands of a finite size are 
thermodynamically stable. Our ah initio results for the 
surface energies and surface stresses, however, indicate 
that neither in the case of InAs/GaAs nor InP/GalnP 
would this mechanism result in the formation of stable 
islands. 

For the most widely studied system of InAs islands 
on GaAs, a range of growth parameters seem to have 
been established for producing islands of certaiii^umfiprm 
densities and a rather narrow size distribution,Ejc3lj al- 
though these growth parameters are still subjects of dis- 
cussion in the literature. Spectroscopic studies of the is- 
lands, after thej« were capped by barrier materials, have 
been reported.cjca To investigate the equilibrium shape 
and stability of the coherent strained InAs islands at low 
island densities, we hav£-£mplied a hybrid method to cal- 
culate the total energy.EETEJ In this approach the energy 
gained by island formation is decomposed in the following 
form: 
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where E'^'^^'^'^ is the gain in deformation energy when the 
material forms a strained island instead of a biaxially 



strained film, iJ^^urface jg ^^^ ^^^^ ^^ surface energy due to 
the creation of facets on the sides of the island instead of 
the surface covered by the base of the island, and iJ^'^s^ 
is the energy cost for the creation of sharp edges. For 
an isolated island to form at all in preference to a film. 
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i? must be negative. The surface energy E^ 
calculated ab initio, applying density-functional theory 
(DFT). In Ref. |o| the surface energies £:«"'-facc corre- 
sponded to those of unstrained surfaces. This approx- 
imation was subsequently improved upon in Ref. |3l|, in 
which the renormalization of the surface energies due to 
surface stress was taken into account. However, the cor- 
rections amounted to a reduction of the surface energies 
by at most 11 % and left the prediction for the equilib- 
rium shape qualitatively unchanged. Furthermore, the 
edge energy E'^'^^'^ was estimated to be negligible, pro- 
vided the island size is not too small. The elastic strain 
field can be treated within a continuum approach. Thus 
^rciax jg calculated within linear elasticity theory using a 
finite-element method (FEM) . The effect of nonlinearity 
was seen to be small. 

As a result, we obtained a volume-dependent optimum 
shape for the InAs islands, which can be described as a 
(OOl)-truncated pyramid with {101}, {111}, and {111} 
faces. However, the diversity of experimentally observed 
island shapes appears not to be reconcilable with the as- 
sumption of thermodynamic equilibrium. Among the ex- 
perimental results, th&,mostly|-sauare- based islands have 
faces {101},Ea {105},EI {113}P {136},Eil and a series of 
islands of low aspect ratios-whose morphologies change 
according to the coverage.E3 We take this difference as 
an indication that a sensitivity to growth conditions and 
kinetic effects have to be featured in a growth theory, in- 
cluding the possibility that the deposited material may 
migrate orr-segregate as witnessed in the ff^Y^ *^^ quan- 
tum wells,E3E3 and self-organized islands.cZTLJ 

While the predictions of the hybrid method for the 
equilibrium shape of InAs/GaAs coherent islands have 
not been borne out so far in experiment, i.e., in the 
'window' of growth condition assumed, the experimen- 
tal characterization of ptte three-dimensional islands in 
the InP/GalnP system,EZl briefly summarized in the next 
paragraph, seems to indicate a better chance for an equi- 
librium approach to be valid. Therefore, InP/GalnP 
represents an excellent benchmark system to show both 
the applicability of the hybrid method and the notion 
that the shape equilibration of coherent islands can be 
achieved. 

The growth of InP/GalnP by metal- 
organic vapor phase epita^f (MOVPE) has consistently 
yielded uncovered islandgl3 that are significantly larger 
(~ 45 X 60 nm^) than their counterparts in InAs/GaAs 
(~ 12 X 12 nm^) grown by molecular beam epitaxy 
(MBE). Transmission electron microscopy (TEM) and 
atomic force microscopy (AFM) have yielded strong 
evideace, that they have facets of low Miller indices 
only,CJ'Ell rather reminiscent of the previous prediction 
for InAs/GaAs. They are stable against annealingtZl of 




several minutes at the growth temperature of 580 °C, and 
their morphology has been reported to repjain unchanged 
after overgrowth with capping materiaLcd Although for 
these materials there is a lack of specific observations of 
segregation inside the islands and diffusion between the 
islands and the barrier, it is probable that both takes 
place. To what extent they affect the shape and size of 
InP/GaP islands remains to be determined. It is also 
becoming clear that the spectroscopy of the islands de- 
pends not only on the volume but also significantly on 
the strain distribution inside the islands, and tiifi-4atter 
is greatly influenced by the shape of the islandsH3a The 
results for InP/GalnP complement those for InAs/GaAs, 
e.g., growth of|-quaiifa|im wells on InP isla 
characteristics,OL3Ei-|photoluminescence,^t 
tical gain and lasingj^ InP islands used 
induce quantum dots in a quantum wellj 
levels formation in InP islands P^.theoretical study of the 
electronic states of the islands.oS 

We are thus encouraged to conduct a study of 
InP/GalnP-islands parallel to our previous work for 
InAs/GaAs, and hence enlarge our general understand- 
ing of coherent strained islands. Furthermore, gaining 
an insight into the strain distribution of capped islands 
would help towards a better understanding of the poten- 
tial experienced by the electrons and holes at the quan- 
tum dot and in its vicinity. Of course, if the effect of 
phase segregation inside the island or diffusion of atoms 
to and from the barrier turns out to be substantial, our 
approach would need to be generalized. In any event, our 
findings could still be used as the starting point for a yet 
more realistic modeling. 

The organization of the work is as follows: First, we 
present the DFT results for the surface energies. In Sec- 
tion III, we describe briefly our FEM simulations and 
derive the equilibrium shape. We shall indicate how Ost- 
wald ripening emerges from our model, assuming that 
no other faster kinetic processes have preempted its time 
development. In Section IV we compare the strain distri- 
butions of uncapped and capped islands. The implication 
for the use of strained islands as a stressor for quantum 
wells will be pointed out. A generalization from InP/GaP 
to InP/GaxIni_xP islands, i.e., to systems with a differ- 
ent lattice mismatch, is included in the final Section V. 



II. SURFACE ENERGY 

As in the case of InAs we assume that the rele- 
vant surface reconstructions, i.e., those of lowest en- 
ergy, correspond to the low-index surface orientations 
{110},{100},{111}, and {111}. The details of the calcu- 
lation of surface energies by DFT are described in Refs. 
^,^. We apply the local density approximation to the 
exchange-correlation energy-functional and use ab initio 
norm-conserving, fully separable pseudopotentials. The 
plane- wave expansion of the wave function has an energy- 



TABLE I. The relaxed surface reconstructions of InP for 
fip — /ip""' — 0.1 eV, and surface energies. 



orientation 


reconstruction 


surface energy 

[meV/A''] 


(110) 


cleavage 


55 


(100) 


/32(2 X 4) 


62 


(111) 


(2 X 2) In-vacancy 


62 


(111) 


(2 X 2) P-trimer 


44 



cutoff of 10 Ry, and the k-summation makes use of a uni- 
form Monkhorst-Pack mesh with a density equivalent to 
64 k-points in the complete (1x1) surface Brillouin zone 
of the (100) surface. We neglect the correction of the sur- 
face energies due to strain, as, for the InAs/GaAs- islands, 
we found it to be small and not affecting mec conclusions 
concerning the island shape and stability.LlI 

The surface atomic structures for different surface ori- 
entations are shown in Fig.n^. The corresponding sur- 
face energies as a function of the phosphorous chemical 
potential are shown in Fig.0, where the left and right 
vertical dashed lines mark the limits for In- and P-rich 
environments, respectively. They are characterized by 
the coexistence of the InP surface with either an In or a 
P bulk phase. Since epitaxial growth takes place mostly 
in a P-rich environment, we list in Table I the surface 
energies of the stable reconstructions for the chemical 
potential ^p = ^p""^ — 0.1 eV. It will be shown in Section 
V that the experimentally observed coherent islands are 
best compared with the theoretical results at this chem- 
ical potential. 

There are a great deal of simLlpdties between InP and 
InAs surfaces in equilibriumE3'E30 and some interesting 
differences. Both the (110) and (111) surface energies 
are independent of the phosphorous chemical potential. 
The (1x1) relaxed cleavage surface and the stoichio- 
metric In vacancy structure are the stable reconstruc- 
tions for the (110) and (111) orientations, respectively. 

The surface energies are 55 and 62 meV/A , respectively. 
These are to be compared with those of InAs of 41 and 

42 meV/A . The cleavage surfaces of the (110) orienta- 
tion in InP have been obseppd experimentally using low- 
energy clectron-diffractioiH (LEED) and for iipth InP 
and InAs using low-energy positron-diffraction.Ej Api)FT 
study of InP(llO) surface has also been carried out.Ej We 
find the equilibrium structure of the (111) surfaces to be 
the same as for InAs. On the (111) surface in P-rich en- 
vironment, the P-trimer reconstruction is preferred, i.e., 
within some interval near the right-hand dashed line in 
Fig.|. At HP = ^^"i'^ - 0.1 eV the InP (111) surface en- 
ergy is 44 meV/ A compared to 36 meV/ A for InAs. A 
(\/T9 X -\/T9) reconstruction of the (111) surfaces seems 
not relevant for our present study, as the area on side 
faces of the island might be too small to accommodate 
this reconstruction without having to invoke edge effects. 
The (100) reconstructed surfaces of InP and InAs had 
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FIG. 1. Atomic structure models for the different InP sur- 
faces, top and side views. Filled and open circles denote In 
and P atoms, respectively. 



been generally assumed to be similar, too. However, re- 
cent DFT studies and experimental observations by sev- 
eral groups using low-energy electron diffraction (LEED), 
reflection anisotropy spectroscopy (RAS), soft x-ray pho- 
toelectron spectroscopy (SXPS), and scanning tunneling 
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FIG. 2. InP surface energies of the (110), (100), (111), and 
(ill) surface orientations as a function of the phosphorous 
chemical potential. The thick lines highlight the calculated 
surface energies of the reconstructions of lowest total energy. 



microscopy (STM) have yielded evidence that there is a 
qualitative difference between the. atomic structm'c of the 
InP and the InAs(lOO) surfaces J£jEZI In a very P-rich en- 



vironment, i.e., at /ip rather close to /^p"''^, DFT predicts 
the c(4 X 4) reconstruction to be more stable than the 
/32(2 X 4) reconstruction, which is not the case for InAs. 
For InP in moderately P-rich and InAs in As-rich envi- 
ronments, /32(2 X 4) is the common stable reconstruction. 
However, in an In-rich environment InP displays a mixed 
dimer reconstruction, in contrast to the a (2 x 4) recon- 
struction for InAs. However, for the growth of strained 
islands in a moderately P-rich atmosphere, it is still the 
/32(2 X 4) reconstruction that is selected as in InAs/GaAs, 

with the surface energy equal to 62 meV/A compared 
to 44 meV/A for InAs. 



III. EQUILIBRIUM ISLAND SHAPES 

For a fixed volume of the island (or, equivalently, a 
fixed number of atoms), the equilibrium shape minimizes 
the total energy of the system, E'^°^'^'^ of Eq. (1), with re- 
spect to all possible shapes. As explained above, the 
surface energy £;>5"''faco jg calculated for the unstrained 
surface and the edge energy _B'='^sc jg neglected.Eil For a 
square-based island bounded solely by the four {101} sur- 
faces, the length of the base is chosen equal to 12.9 nm 
and the height to 6.5 nm. This determines the volume 
that remains unchanged for all calculated islands. The 
island dimensions vary only moderately for truncated is- 
lands bounded by the other low-index surface planes we 
have included in this study. By construction the islands 
are placed 40 nm apart, and reside on a substrate of thick- 
ness 24 nm. As in the case of InAs,EEll£^ the thickness of 
the wetting layer has been set to zero. ■— ■ 

We have used the commercial product MARCE3 to 
perform the finite-element simulations. One notable fea- 
ture of all commercial products is that the preferred fi- 
nite element (FE) is an eight-node hexahedron or three- 
dimensional arbitrarily distorted cube. It has been shown 
that it is superior to the simple tetrahedron in terms of 
fast convergence and computing speed. We have also 
adopted this class of FE for all our calculations, with the 
understanding that corners of the hexahedron can merge 
to form wedge-shaped or pyramid-shaped elements. 

Throughout this work we have taken InP to be the 
island and wetting layer material, and GaP to be the 
barrier material, or the substrate when the island is un- 
capped. The correspondJBg experimental lattice con- 
stants and elastic moduliEj are listed in Table II. The 
lattice mismatch is a :— (acaP — ainp)/aniP = —7.1%. 
Although this mismatch is quite large, we restrict our- 



TABLE II. The experimental lattice constants 
tic moduh cn, Ci2 and C44 of InP and GaP. 


a and elas- 


a A cii [GPa] C12 [GPa] 


C44 [GPa] 


InP 5.87 102 58 
GaP 5.45 141 63 


46 
72 




FIG. 3. The relaxation energy per unit volume E''°'^'^^ /V 
versus £;su''face^^2/3 £qj. YnP islands. Square : square-based 
island with four {101} facets. Diamond : square-based is- 
land with two {111} and two {III} facets. Triangles up 
huts with two {111} and two {III} facets. Triangles down 
square-based {101} island with {111} truncated edges. Dots 
islands with four {101}, two {111}, and two {111} facets. The 
dashed hne is the hne of constant total energy ^j^eiax _|_ ^surface 
= constant that selects the equilibrium shape for the volume 
V = 3 X lO^A . See text for further explanations. 



selves to linear elasticity theory. Taking the elastic mod- 
uh from Table II, the deformation energy per unit volume 
of the biaxially strained uniform InP film amounts to 3.0 
meV/Al 

In our FEM calculations of E'^'^^^^ the number of FE 
varies slightly, depending on the island's bounding sur- 
faces. The number of FE is increased in the areas 
where the deformation energy density is large until a 
1 % accuracy is achieved. In general, there are approxi- 
mately 7000 FE in total, of which approximately 2000 are 
distributed in the island. For volume-conserving trun- 
cated islands, the elastic energy can be accurately ap- 
proxima±£|d from the untruncated value by an analytical 
formula.Lj 

The anisotropy of the surface energy driving the for- 
mation of particular crystal facets is an essential aspect 
of our approachJiJ We have investigated island shapes 
which are bounded by the low Miller-index surfaces de- 
scribed in Section II. The collection of shapes we have 
conside|C£|d_-.is the same as in the study of the InAs 
islandsEj'Ell The island base has the orientation (001) 
and the top of the island may be truncated by a plane of 
the same orientation (see Fig. 6 of Ref. pl|) . 

The results of energy evaluated according En. (1) are 
summarized in a scale- invariant manner in Fig.^, where 
the ordinate x = £;surtaco^y2/3 ^^^^ abscissa y — 
j^rcinx^Y^ We have taken the surface energies at the 



phosphorous chemical potential /xp = /Up"''^ — 0.1 eV. 
The solid symbols denote untruncated islands for which 
a full FE calculation has been carried out. Qualitatively 
these results are similar to those of InAs/GaAs.E]'Eil We 
find that a square-based island with {lOlj-facets only 
(solid square) has a larger bulk deformation energy than 
a square-based island with two {111}- and two {111}- 
facets (rhombus). The latter has steeper side facets which 
allow a more efficient stress relaxation. The line that em- 
anates from each of the solid symbols joins up the small 
dots for which E'^^^'^^ was derived from the analytical ap- 
proximation for the volume-conserving truncated island. 
For a given volume V, islands with the same total energy 
^totai j-g Qj^ ^Yie straight hne 
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< j^Y /-^ . X — constant, 
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plotted as a dashed line in Fig.0. (Note that the edge en- 
ergies iJ^'^sc have been neglected.) It is clear from Eq. (g) 
that the volume of the island V is related to the negative 
slope of the line. For a given volume V, the equilibrium 
shape of the island is determined by the first point of 
contact from below of the straight line, Eq. (||), with the 
calculated island-energy curves. 

For a given island shape, E'^'^^^^ and iJ^^urfacc g^g^jg jjj^g 
V and y^/"^, respectively. This means that the data in 
Fig. are invariant against a change of the volume of 
all islands by the same factor. However, the slope of the 
straight line, Eq. (||), changes, and hence a different equi- 
librium shape will be selected. By inspection, the results 
of Fig.0 indicate that for larger V, i.e., smaller negative 
slope, untruncated islands are most likely to be selected, 
while smaller V favors truncated islands as equilibrium 
shape. The interplay of i^'^'^'^ and £;s">-face j^ Eq. (1), and 
how their relative weightings change as a function of the 
island volume have been discussed and supported by ex- 
plicit calculations in Ref. |3y. Fig.y serves as a compact 
way of encapsulating the result for the equilibrium shapes 
for arbitrary volumes, provided that the volume (or the 
number of atoms) of the island is not too small to inval- 
idate our hybrid ansatz, Eq. (1), or too large such that 
the island is no longer coherent (dislocations appear). 

We display in Fig.H examples of equilibrium shapes 
for three different island volumes. They illustrate nicely 
that for a smaller volume a relatively larger amount of 
material is truncated from the top of the island. Ad- 
ditionally, we can also inspect the difference in shapes 
due to a variation of the chemical potential: The {111}- 
facet, which is favored in a P-rich environment, is con- 
sistently more prominent for /ip = /ip""^ — 0.1 eV than 
for /ip = /Zp""^ — 0.2 eV. The equilibrium shapes shown 
in Fig. y have m uch aimilarities with the experimentally 
observed islands ,c3cil in particular all observed facets are 
well accounted for by our theory. Note that the exper- 
iments have been carried out for InP/GalnP instead of 
InP/GaP, thus the experimental and theoretical lattice 
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FIG. 4. The equilibrium shapes of coherent strained InP 
islands for three different volumes and two different chemical 
potentials. The unstrained volumes in ascending order are 
approximately 9 x fO''A^ 3 x 10^ A^ 6 x fO^A^ 



mismatch differ by nearly a factor of two. However, as 
will be discussed below, this can be accounted for by a 
simple re-scaling of the island volunQ£s. 

In the theory of Shchukin et ai,E3 a combination of 
material parameters determines whether thermodynamic 
equilibration drives the system towards a stable array of 
islands or whether Ostwald ripening occurs. A stable 
array would arise, if E^^'^^'^'^ in Eq. (1) is modified by the 
surface strain in such a way that it changes sign. Making 
use of our quantitative surface energies, we can estimate, 
whether this condition is likely to be fulfilled in case of 
InP or not. 

The change in surface energy i<;>5"'facc jg given by a sum 
over all facets of the island: 
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^base g^^^ ^i denote the areas of the base and the i- 
th facet on the surface of the island, respectively, 9^ is 
the angle between the i-th facet and the substrate sur- 
face, and v"^ is the ratio of A^ projected onto the base 
and A!^^"^ , The 7' are the surface energies from Table 
I. Each term within the summation of Eq. (|^) is posi- 
tive, and hence E^™^'^^'^ is positive. This has the con- 
sequence that, for a fixed equilibrium island shape, the 
volume derivative of Eq. (||) is always negative, i.e., the 
total energy per unit volume decreases as the volume in- 
creases, thus favoring larger and larger islands, leading 
to Ostwald ripening. To change the sign of the respec- 
tive individual contribution to Eq. (0) , a strain-induced 
renormalization of the surface energy of approximately 



20 % for the facets with (101) and (111) orientations, and 
40% for the facet with (111) orientation would be nec- 
essary. For the case of InAsyCaAs, the renormalization 
amounted to at most 11 %.Eil Therefore, it seems reason- 
able to argue that Ostwald ripening is also the fate of the 
InP islands. Furthermore, a majority of the InP islands 
are grown on GalnP-substrates whose lattice mismatch 
is about 50 % smaller than in the InP/GaP system. This 
further diminishes the importance of the renormalization 
of the surface energy due to the strain field. Altogether 
this discussion seems to indicate that we are outside the 
parameter range for which arrays of stable islands occur. 
It is worth noting that— the parameter range that con- 
forms to Shchukin et aJ,lla implies that £;«">-faco j^ Eq. (1) 
changes sign due to stress, i.e., one has the unusual situ- 
ation that there is a gain in total energy by creating new 
surfaces. Furthermore, Shchukins model includes edge 
energies and island-island interaction to prevent the is= 
lands from coalescing. In comparison, a new mcthodcj 
to tackle this problem of the stability of the islands has 
a much less restrictive ansatz and the agreement with 
experiments is encouraging. 



IV. COMPARISON OF UNCAPPED AND 
CAPPED ISLANDS 

In this section we compare the strain fields for capped 
and uncapped islands. It is well documented that during 
overgrowth of GaAs over InAs islands there is signifi- 
cant migration of atona^-as well as change in the mor- 
phology of the islands, LjlJ and that In segregates during 
the OYeffirpadiLjaf GaAs over InGaAs islands or quantum 
wells. OEElH'HiZfl A well-studied sample of InP strained 
islands was,«jown on InGaP lattice-matched to a GaAs 
substrate.E^ES These InP islands seem quite stable - when 
uncapped, they survive a few minutes of annealing at 580 
°C without changes and retain their shapes after over- 
growth. These may be signs that the MOVPE-grown InP 
islands are less affected by atom migration and In segre- 
gation during overgrowth than the MBE-grown InAs is- 
lands. Therefore, until more specific observations of these 
effects are available, it would be interesting to compare 
the strain distributions of uncapped and capped islands 
without taking into account any change of shape or In 
segregation. While we are aware that the present results 
may be subjected to revision, we shall discuss below some 
rather general effect of the capping layer on the strain 
distributions. 

We have adopted the shape and size o£-the island that 
had been observed by Georgsson et ai,E3 displayed in 
their Fig. 1, i.e., a truncated pyramid grown on the (001) 
plane, with the defining planes for the side walls being 
(Oil), (101), (111), (Oil), (101), and (111). The base is 
elongated in the [110] direction and has a length of 66 nm, 
a width of 40 nm in the [110] direction, and the height 
is 18.4 nm. The island resides on a 0.6 nm thick wetting 
layer. By construction the nearest-neighbor islands are 
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FIG. 5. Strain distribution of an uncapped truncated InP 
island on a wetting layer and a GaP substrate. 
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140 nm apart along the [110] and [110] directions. The 
size of our ceh is chosen large enough such that away from 
the island and the wetting layer the barrier material is 
not significantly strained. Since the island has reflection 
symmetry across two planes, only one quarter of the unit 
cell needs to be explicitly considered in the simulation. 
In the case of the uncapped island we have used in total 
1872 FE, while for the capped island we used 3756 FE. 
Our simulation of InP/GaP is not exactly a reproduction 
of the InP/GalnP coherent island in Ref. ^ because the 
lattice mismatch of the latter was smaller, a « 3.7 %. We 
discuss the scaling with respect to the lattice mismatch 
at the end of this section. 

We display in Fig.|^ two views of the trace of the strain 
tensor field, en + 622 + 633, for an uncapped island. In 
Fig. H (a) a part of the wetting layer has been peeled off to 
reveal the strain distribution in the substrate. Qualita- 
tively the results are quite skpilar to those for InAs/GaAs 
without the wetting layer.Ell One notes that the unhin- 
dered relaxation of the island in the upward direction has 
produced a small area of tensile strain at the top of the 
island. This seems to be a common result in all our cal- 
culations, whether the island is truncated or not. In the 
substrate, the strain distribution is characterized by com- 
pressive strain around the edges of the island, and tensile 
strain directly underneath the island. Such contrasting 
strains created by coherent strained islands of InP have 
been utilized as strepspp-|to achieve three-dimensional 
quantum confinementE3cZI in sub-surface quantum-well 
structures, and a theoretical strain analysis similar to 
ours by FEM has been reported. In this discussion, 
an InGaAs quantum well (QW) is sandwiched between 
GaAs, with self-organized InP islands grown on top. The 
strain distribution thus produced modulates the conduc- 
tion band of the QW, and creates a potential well. Pho- 
toluminescence spectra of this type of potential well have 
been measured and analyzed, hence this object has also 



FIG. 6. Strain distribution of a truncated InP island on a 
wetting layer, embedded in GaP. 



earned the label of quantum dot. 

For comparison we display in Fig.o two views of a 
capped InP island. In Fig.g (a) part of the capping mate- 
rial is removed to reveal the island and the wetting layer. 
Part of the latter is also peeled off to show the substrate 
underneath. The island-|in this instance is entirely un- 
der compressive strain.O The area of tensile strain pre- 
viously found in the substrate directly underneath the 
uncapped island has been replaced by an area of com- 
pressive strain. Small pockets of tensile strain are to be 
found in the vicinities of sharp corners in the capping 
layer. The disappearance of the pocket of tensile strain 
in the substrate has great significance for the creation of 
quantum dots in QWs as described above. The amount 
of capping material affects the strain-induced potential 
well in the QW. This phenomenon is quite general and 
not limited to strained InP islands. If we assume that the 
strain interaction is given by a 'hydrostatic' potential, 
i.e., V oc tre, our present calculation shows the extreme 
case of destroying this potential well completely. 

A closer examination of the components en, €22, 
and £33 directly underneath the island reveals that en 
and 622, are tensile and remain qualitatively unchanged 
whether the island is capped or not. On the other hand, 
£33 is slightly tensile for the uncapped island and over- 
whelmingly compressive for the capped island. These me- 
chanical responses seem understandable. For capped or 
uncapped islands the cause of the biaxial tensile strain in 
the substrate perpendicular to the [001] direction is the 
same. However, the uncapped island heaves up in the 
[001] direction because the relaxation in this direction is 
unhindered. For the capped island, the capping material 
prevents the unrestrained heaving in the [001] direction. 
Hence the island exerts a considerable pressure onto the 



substrate as well, i.e., material in the substrate is com- 
pressed. Therefore the thickness of the capping layer has 
a great influence on the strain in the substrate. 



V. DISCUSSION 

The experimental data on the growth of InP islands 
mostly involve I nB grown , on GalnP lattice- matched to a 
GaAs substrate.E^lci'cSrEil In the cases where InP islands 
were used as stressors of QWs, they were grown on GaAs. 
In all these cases the lattice mismatch amounts to about 
-3.7%. TEM and AFM have been used to characterize 
the coherent strained InP islands and their lateral sizes 
are usually found to be roughly 40 nm x 60 nm and the 
height between 15 nm to 20 nm. 

In our study we have assumed a GaP substrate with a 
lattice mismatch of —7.1 %. Therefore, our data have to 
be rescaled before comparing our equilibrium shapes to 
the experimental observations. The equilibrium shape at 
a given volume V and lattice mismatch a is calculated by 
minimizing the total energy i?'°*'''(a, V,/) with respect 
to the island shape /. As above, we omit the contribu- 
tion from the edge energy and the renormalization of the 
surface energy due to the strain, and make use of the 
scaling behavior of the elastic energy E'^°^^^ and surface 
energy £;surtaco^ 






relax 



■surface 






(4) 
(5) 



where Eq'^^^^{I) and Eq^'^^'^'^^{I) are constants of propor- 
tionality that depend only on the shape. The total energy 
of the equilibrium island follows from 



E 



total 



(a,y) = 



a t^min 



^rolax^j) -H £;j5"f^'='=(/) a~ V~l/3 



(6) 



Therefore the optimum island shape is the same for vol- 
umes V,V' and lattice mismatches a, a', respectively, 
provided that the relation 



'-V^/^ = a'V^/^ 



(7) 



holds. This simple relation will be altered, however, if the 
different elastic moduli for GaP and GalnP are taken into 
account. 

Let us select from Fig.^ the theoretical equilibrium 
island shape pertaining to the chemical potential fip = 
^buik _ 0.1 eV, and volume V^ « 3 x 10^ A , because it 
bears the closest resemblance in shape to the island ex- 
amined by AFM in Ref. 41, We assume the volume of 



the latter to have the dimension given in Section IV, i.e., 

V = 3.7x 10^ A . Therefore, for an island of volume 

V ^ 3 X lO^A , formed with a mismatch of —7.1 %, to 

transform to an island of volume V = 3.7 x 10^ A , the 
mismatch, as given by Eq. (Q), ought to be —4.6%. In 



view of all the uncertainties this appears to be reason- 
able. Note that the interdiffusion of In and Ga from the 
InP island and GalnP substrate, would decrease the as- 
sumed lattice mismatch of —3.7%, hence worsening the 
agreement with our estimate, and would also alter the 
effective elastic moduli of the island and barrier materi- 
als. Although this topic has yetJjO be explored in these 
materials, there is clear evidencdZj for In and Ga interdif- 
fusion in the growth of InxGai_xAs islands embedded in 
GaAs. A quantitative clarification of interdiffusion seems 
to be crucial for a better understanding of the equilibrium 
shape and stability of the strained islands. In principle, 
a concentration profile or phase separation in the island 
can be incorporated into the FEM calculation. 

In conclusion, we have reported an application of the 
hybrid method to study the equilibrium shape of coherent 
strained islands of InP on GaP or GalnP substrates. Our 
results are in qualitative agreement with experiments, in 
that all the observed facets are accounted for by our the- 
ory. We take this as an encouraging evidence that the 
coherent InP islands most often observed in experiments 
are close to local thermodynamic equilibrium. This is in 
contrasts to the apparent discrepancy between the cal- 
culated and observed island shapes for the InAs/GaAs 
islands. In view of the distinct variety of observed InAs 
island shapes this discussion points towards the impor- 
tance of kinetic effects as the missing ingredient of a more 
comprehensive theory for those islands. The contribution 
of the surface energies to the total energy of an InP is- 
land is positive, and it seems unlikely that the sign of 
this contribution would change, if effects due to surface 
stress and strain were included. Thus our results favor 
the Ostwald ripening of the islands (as opposed to the 
creation of a thermodynamically stable array of islands). 
Finally, we have compared the strain fields of uncapped 
and capped islands. The capping material distinctly af- 
fects the strain fields in the substrate. This is of relevance 
for quantum dots induced in quantum wells by stressors. 
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